Final Results from the KTeV Experiment on the Decay K h — ► 7r°77 
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We report on a new measurement of the branching ratio TS(Kl — > 7r°77) using the KTeV detector. 
We reconstruct 1982 events with an estimated background of 608, that results in B(Kl — > 7r°77) 
= (1.29 ± 0.03 s tat ± 0.05 S yst) x 10" 6 . We also measure the parameter, av, which characterizes the 
strength of vector meson exchange terms in this decay. We find av = —0.31 ± 0.05 s tat ± 0.07 S i, S i. 
These results utilize the full KTeV data set collected from 1997 to 2000 and supersede earlier KTeV 
measurements of the branching ratio and av- 



PACS numbers: 13.20.Eb, 11.30.Er, 12.39.Fe, 13.40.Gp 
I. INTRODUCTION 

The decay Kl — ► 7r°77 provides important checks of 
low-energy theories of strange meson decays. In Chiral 
Perturbation Theory (ChPT) the branching ratio for this 
decay can be determined with no free parameters up to 
0(p i ). However, the first measurements of the branch- 
ing ratio for Kl — > 7r°77 [l|, Q, Q were approximately 
three times larger than the predicted 0(p 4 ) branching 
ratio of 0.68 x 10 _6 [J|. Extending the theory to 0(p 6 ) 
and including vector meson exchange terms raise the 
branching ratioprediction to be consistent with the mea- 
sured values [HQ. The vector meson contributions can be 
parametrized by an effective coupling constant ay- Non 
perturbative calculations for the Kr, — > 7r°77 rate have 
also been performed^- 

The Kl — > ir°jj decay is important also because of 
its implications for the related decay Kl —* ir°l + l~, 
where l + l~ can be either e + e _ or Currently, 
the best limits for these decays are B(Kl — > w°e + e~) 
< 2.8 x 10- 10 [| and B(K L -» ttVaO < 3.8 x 10" 10 i, 
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both at the 90% confidence level. The ex pec ted branch- 
ing ratios are approximately 1 — 43 x 10 11 [lot [TT|| . There 
are three contributions to the Kl — > tt°1 + 1~ decay, clas- 
sified in terms of their CP symmetry; one conserves CP 
symmetry, one violates it indirectly and one directly. The 
direct CP violating amplitude is of interest within the 
Standard Model but also can show signs of new physics 
pl| . leading to an enhancement of the Kl — > rate. 
In order to determine the direct CP violating terms, one 
must first determine the other two amplitudes. The indi- 
rect CP violating amplitude can be determined from the 
decay K$ — ► ir°l + l~ , and the NA48 experiment has mea- 
sured B(K S -> 7r°e+e-) = 5.8±|:| x lO" 9 ^. Because 
the Kl — ► tt°I + I~ decay can proceed via a CP conserv- 
ing two-photon exchange, the CP conserving terms can 
be probed using Kl —> tt°jj. A precise measurement of 
the parameter ay can be used to determine the size of 
the CP conserving amplitude in Kl — > ir°l + l~. 

There have been a number of previous measurements 
of K L -> 7r°77 from the E731, NA31, NA48 and KTeV 
experiments. [H, S, 1 El El- 

The two most recent 
measurements are the NA48 result of (1.36 ± 0.03 s t a t ± 
0.03 syst ± 0.03„ orm ) x 10" 6 and the KTeV result of 
(1.68 ± 0.07 stat ± 0.08 syst ) x 10~ 6 . Both of these results 
are significantly more precise than the E731 and NA31 
results. However, the NA48 and KTeV results differ by 
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nearly three standard deviations. The measurement dis- 
cussed here supersedes the previous KTeV result and rec- 
onciles the difference between these two results. 



FIG. 2: Transverse view of the KTeV Csl calorimeter. The 
smaller blocks are located in the central region with the larger 
blocks located in the outer region. 



II. THE KTEV DETECTOR 

Data used in this analysis were collected during three 
running periods in 1996, 1997 and 1999 using the KTeV 
detector at Fermilab. Because of the similar topology 
between the Kl — ► 7r°77 d ecay s and the Kl — > 7r°7r° 
decays used to measure e'/e [HI], we recorded the Kl — > 
7r°77 events during the same collection period used for 
the KTeV e'/e measurement. 

The KTeV experiment[l6] is a fixed-target experiment 
built to study decays of neutral kaons. A schematic of 
the detector is shown in Figure [TJ Two neutral kaon 
beams were produced through interactions of 800 GeV/c 
protons in a 30 cm long beryllium oxide target. The 
resulting neutral particles passed through a series of col- 
limators and absorbers to produce two nearly parallel 
beams. Charged particles were removed from the beams 
by sweeping magnets located downstream of the collima- 
tors. A vacuum decay volume extended from 94 to 159 
meters downstream of the target, and was far enough 
away from the target that the vast majority of the K$ 
component had decayed away. An active regenerator was 
located within the vacuum region, approximately 123 
meters downstream of the target. This regenerator al- 
ternated between the two neutral beams to generate a 
Kg component in one of the beams. The beam that co- 
incided with the regenerator was called the regenerator 
beam, while the other beam was denoted the vacuum 
beam. For this analysis, we only considered decays from 
the vacuum beam. To reject photons, primarily from de- 
cays of Kl — > 7r 7r°7T , the decay volume was surrounded 
by photon veto detectors, that rejected photons produced 
at angles greater than 100 milliradians with laboratory 
energies greater than 100 MeV. A kevlar and mylar vac- 
uum window with a radiation length of 0.14% covered 
the downstream end of the vacuum decay region. 

The most critical detector element in this analysis was 
the pure Csl electromagnetic calorimeter |15j. The Csl 



calorimeter, shown in Figure [H was composed of 3100 
blocks in a 1.9 m by 1.9 m array with a depth of 50 cm 
corresponding to 27 radiation lengths. Two 15 cm by 15 
cm holes were located near the center of the array for the 
passage of the neutral beams. For photons with energies 
between 2 and 60 GeV, the calorimeter energy resolution 
was below 1% and the nonlinearity was less than 0.5% 
per 100 GeV. The position resolution of the calorimeter 
was approximately 1 mm. 

Two levels of hardware triggers were used in the KTeV 
experiment. For the Kl — > tt°jj events, the first level 
trigger required the event to deposit more than approxi- 
mately 25 GeV in the Csl calorimeter with less than 100 
MeV in any of the photon vetoes. The second level trig- 
ger utilized a hardware cluster processor that counted the 
number of separate clusters in the Csl calorimeter [l7| . 
Each cluster had to have an energy greater than approx- 
imately 1.0 GeV and the total number of clusters in the 
Csl calorimeter was required to be equal to four. 

Events that satisfied the hardware triggers were also 
required to satisfy an online software filter. This filter 
required that one of the six possible combinations of two 
photons reconstructed near the ir° mass. In addition, 
the filter required that the decay vertex reconstructed 
upstream of 155 m for the 1996-1997 data and 140 m for 
the 1999 data. This requirement was tightened for the 
1999 data to help reduce the trigger rate for the Kl, — > 
7r°77 sample. These trigger requirements also selected 
Kl — » 7r°7r° events which were used as a normalization 
mode for calculating the Kl — > 7r°77 branching ratio. 

III. EVENT RECONSTRUCTION 

The Kl — > 7r°77 final state consists of four photons 
with no other activity in the detector. In this analysis 
we require that all events have exactly four clusters in 
the Csl calorimeter and that the energy of each cluster 
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is greater than 2.0 GeV. To reduce contamination from 
events originating from the regenerator beam, the center 
of energy (Eq. 9 of [l5.]) is required to be within the Csl 
beam hole corresponding to the vacuum beam. 

In the decay Kl — > 7r°77 the positions and energies of 
the four photons do not provide enough constraints to de- 
termine both the decay position and the invariant mass of 
the system. Therefore, we assume that the four-photon 
invariant mass is equal to the kaon mass, and reconstruct 
the decay vertex position (z) from the calorimeter infor- 
mation. For a 7T° decaying into two photons, one can 
determine the two-photon mass using the following rela- 
tion: 



m 12 



(1) 



where E\ and E2 are the energies of the two photons, 
ri2 is the distance between the two photons at the Csl 
calorimeter, and Azcsi is the distance between the decay 
vertex and the Csl calorimeter. Using the position of the 
reconstructed decay vertex, we determine the two-photon 
mass for each of the six possible combinations and choose 
the combination with the reconstructed mass closest to 
the known ir° mass. If the closest mass combination dif- 
fers by more than 3 MeV/c 2 from the known ir° mass, 
we reject the event. The total energy of the kaon sys- 
tem, determined from summing the energies of the four 
clusters, is required to be between 40 and 160 GeV. Af- 
ter these requirements the data sample is dominated by 
backgrounds from Kl — > 7r°7r°7r and Kl — > 7r°7r° de- 
cays. Additional cuts described in Section [V] are used to 
reduce these backgrounds. 



IV. MONTE CARLO SIMULATION 

A detailed Monte Carlo simulation was used to esti- 
mate the detector acceptance and the background level 
in our final sample. Our Monte Carlo simulates the kaon 
production at the target and propagates the kaon am- 
plitude through the detector. The kaon then decays ac- 
cording to the appropriate decay mode, and the resulting 
daughter particles are traced through the KTeV detector. 
The interaction of the decay products with the detector 
is simulated and the detector response is then digitized. 

Details of the simulation for all detector components 
are given in [l5| ; here we focus on the simulation of the 
Csl calorimeter. To simulate the response of the Csl 
calorimeter, we used a library of photons generated us- 
ing GEANT simulations [18]. The library contained in- 
formation deposited into a 13 x 13 array of Csl crystals. 
The wrapping and shims separating each crystal was in- 
cluded in these simulations. This library was binned as a 
function of the energy and position of the incident pho- 
ton. We stored the energy depositions for each crystal 
in 10 longitudinal bins to include the effects of nonlinear 
response along the length of the crystal. 



During the course of our studies we found that the 
GEANT-based shower library was not adequate for de- 
scribing the transverse distribution of the energy in a 
electromagnetic shower. As noted below, we make use 
of this transverse shape to help reduce backgrounds from 
Kl — > 7r 7r°7r° decays. To better simulate the shower 
shapes, we also implemented a data-based shower library. 
These showers were extracted from Kl — > ir°ir° events 
taken during special, low-intensity runs to reduce the ef- 
fects of accidental activity in the Csl calorimeter. The 
data-based shower library was also binned as a function 
of the incident photon energy and the incident position. 

To characterize the transverse energy deposition of a 
electromagnetic shower, we devised a photon shape \ 2 
variable. This variable compares the energy in the cen- 
tral 3x3 crystals of a cluster to the expected energy 
distribution. While a 7 x 7 array of crystals is used for 
accurate cluster energy reconstruction, the photon shape 
\ 2 is determined from the central 3x3 crystals to min- 
imize any biases from accidental activity. As shown in 
Figure [31 Monte Carlo events utilizing the data-based 
shower library match the data better compared to the 
GEANT-based shower library. For our Monte Carlo sam- 
ples, we utilized both shower libraries. The GEANT- 
based shower library was used to determine the energy 
and position of the cluster, while the data-based shower 
library was used for extracting the transverse shape in- 
formation. 



BACKGROUNDS TO K L 



7T°77 



After the event reconstruction discussed in Section [1111 
large backgrounds remain in our data sample. Here we 
discuss the additional criteria used to reduce these back- 
grounds. The major backgrounds in our data sample 
consist of events with neutral beam particles interacting 
in the vacuum window, kaon decays with charged tracks, 
Kl —* 7r°7r° decays, and Kl — > 7r°7r 7r decays, with 
the Kl — ► 7r 7r°7r° decays being the most difficult to re- 
move. Vacuum window interactions can produce 7r°7r° 
and 7r°77 pairs. To remove the vacuum window interac- 
tions, we loop over the six possible two photon combina- 
tions and determine the two-photon decay vertex assum- 
ing the photons resulted from a ir° decay. For each of 
the six possible combinations, we reject the event if the 
decay vertex is downstream of the vacuum window and 
the invariant mass of the other 77 combination is near 
the neutral pion or 77 mass. Events with charged tracks 
are removed by requiring that the total number of hits 
in the drift chamber system is less than 24; a two-track 
event will produce 32 hits in the drift chambers. 

The Kl — * ir°ir° events are easily identifiable because 
both 77 pairs will reconstruct with m 
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, w o mass. Al- 
most all of these events are removed by rejecting events 
in which both 77 masses (mi2 and TO34) are near the 
7T° mass. 77112 is the two-photon invariant mass clos- 
est to the 7T° mass, while TO34 is the invariant mass of 



4 





x 10 


CM 


1200 


O 




i§ 

c 


1000 


CD 




> 
LU 


800 




600 




400 




200 









x 10 


CM 


1200 


O 




c/i 
c: 


1000 







> 
LU 


800 




600 




400 




200 








a ) • 27t°Data 

2ji° MC (GEANT-based) 






' b) 






• 2ti Data 






— 2ti MC (data-based) 


• , J, , ! , 




— — ■ . . , . i . , . , i , , , , i , , , , 



1 2 3 4 5 6 7 8 

Shape x 2 



3000 




All 4y Candidates 

27i° and 7c yy Candidates 




2500 


OA 




Cut f 


2000 








1500 








1000 








500 








ft 









110 112.5 115 117.5 120 122.5 125 127.5 130 132.5 
Decay Vertex m 



FIG. 3: The photon shape \ 2 variable, a) The photon shape 
X variable for K L — > with the dots representing the 

data and the solid histogram the Monte Carlo simulation. The 
Monte Carlo simulation was generated using our GEANT- 
based shower library, b) The photon shape x 2 variable for 
Kl — > 7r°7r events. The dots are the data and the histogram 
represents the Monte Carlo simulation using our data-based 
shower library. The data-based Monte Carlo simulation shows 
marked improvement over the GEANT-based shower library. 



the other pair of photons. In about two percent of the 
Kl — > 7r°7r° events, our choice for 77712 and 77134 did not 
correctly choose both tt° — > 77 decays, and so the cut 
to remove the Kl — > 7r°7r° background fails. To remove 
these events, we also examine the other two possible com- 
binations of the four photons and discard any event in 
which both the 77712 and 77734 values are near the mass of 
the 7T°. 

Because we required exactly four photons, Kl — > 
7r 7r°7r decays can only contribute to the background 
if some of the photons miss the calorimeter or two or 
more photons "fuse" together in the calorimeter. To re- 
duce backgrounds from decays with missing photons, we 
remove events with any significant energy in any of the 
photon vetoes. Also, by restricting the decay region to 
115 < z < 128, we reduce the Kl — > 7r 7r°7r background 
significantly because events with missing photons tend 
to have a reconstructed decay vertex downstream of the 
true decay position. As shown in FigureJH the decay ver- 
tex distribution for Kl — > 7r°7r° events is relatively flat 
downstream of 120 meters. However, the Kl — > 7r 7r 7r° 
background rises sharply as the decay vertex position in- 
creases. 

After applying the cuts described above, there still re- 
mains a significant number of Kl — > 7r°7r 7r° decays; far 
more than the signal from Kl — > tt°jj. These decays 



FIG. 4: The reconstructed decay vertex position from four 
photons for all events (histogram) prior to imposing the decay 
vertex cut. All cuts except for the decay vertex and photon 
shape \ 2 have been applied. The shaded histogram indicates 
Kl — ► 7r°7r° and Kl — ► tt°77 candidates. The rise at large 
values is due to Kl — > 7r°7r 7r decays in which one or more 
photons misses the Csl calorimeter. The position of the cut 
is indicated by the two vertical lines. 

result primarily from events in which two of the Csl clus- 
ters come from fused photons. To remove these events, 
we select events with a small photon shape % 2 . For non- 
fused clusters this variable peaks near zero, while for 
fused clusters this shape x 2 variable becomes quite large. 
Figure [5] shows this variable for both the data and for 
the Kl — * 7r°7r 7r° background. We require the shape x 2 
to be less than 1.8. This cut was chosen to maximize the 
signal significance. For Kl — * 7r 7r 7r background events 
we verified the photon shape x 2 distribution in the signal 
region by reweighting events from the tails of the ir° mass 
distribution and found the resulting shape to correspond 
to our Monte Carlo prediction. The photon shape x 2 is 
our final cut, and reduces our background to a reasonable 
level. 



VI. BRANCHING RATIO AND a v 
DETERMINATION 

After applying all of the cuts described above we find 
1982 events before subtracting background. The final 
777,34 mass distribution is shown in Figure[6l with the data 
well-represented by the signal plus background Monte 
Carlo simulation. The background comprises approxi- 
mately 30% of the total event sample. 

To determine the Kl — > 7r°77 branching fraction, we 
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FIG. 5: The photon shape x 2 variable. The solid histogram 
shows the Kl — > 7r°7r 7r° Monte Carlo while the dotted his- 
togram is the data. The Kl — * 7r°77 signal Monte Carlo 
simulation is the shaded histogram. The sum of the signal 
plus background Monte Carlo simulation agrees well with the 
data. 
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FIG. 6: The 771,34 mass distribution for after all selection cuts. 
The dots show the data, the histogram the sum of the sig- 
nal and background Monte Carlo samples, and the shaded 
histogram the scaled background Monte Claro simulation. 



FIG. 7: The a) energy and b) decay vertex for Kl — > 7r°7r° 
decays. The dots are the data and the histogram is the Monte 
Carlo simulation. 



use the following expression 

B = ((Ntot -N bkg )/N 2 ^ ) x (e 27r o/e w o 77 ) 

x B(K L -> ttV) x B(tt° -> 77), (2) 

where N to t is the number of candidate events, Nbkg is 
the number of background events, A^ 27r o is the number 
of normalization events, and e 27r o and e„.o 77 are the ac- 
ceptances of the Kl — > 7r°7r° and Kl — > 7r°77 events, 
respectively. The acceptances were determined using 
our Monte Carlo simulation described above. The value 



B(K L 



) is the measured Kl — > branch- 



ing ratio. In the previous KTeV analysis, the value of 
B(K L -» 7r°7r°) used was (9.36 ±0.2) x 10~ 4 . We are now 
using the most recent determination of B(Kl — > 7r°7r°) 
= (8.69 ± 0.08) x 10 4 [H, To determine the num- 
ber of Kl — * 7r°7r° normalization decays, we count the 
number of events in the region, 0.130 GeV/c 2 < 77134 < 
0.140 GeV/c 2 . The kaon energy and decay vertex for 
our normalization mode are shown in Figure [7] There 
is good agreement between the data and Monte Carlo 
simulation. 

The numbers used for the branching ratio determina- 
tion are shown in Table |U Note that the acceptances for 
the signal and normalization modes are nearly identical; 
this helps to significantly reduce the systematic uncer- 
tainties due to the acceptance calculation. 

We also extract the value of ay, using the model 
described in [5], from our data by performing a 
two-dimensional maximum likelihood fit to the two 
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Parameter 


1996-1997 


1999 


Total 


Ntot 


989 


993 


1982 


Ntot - N bkg 


670.6 


703.8 


1374.4 


N 2v o Events 


482027 437305 


919332 


Signal Acceptance 


0.0330 


0.0261 


0.030 


Norm Acceptance 


0.0328 


0.0257 


0.030 


K L -> 3tt° Bkg 


313 


288 


601 


K L -» 2tt° Bkg 


5.4 


1.2 


6.6 



TABLE I: Values used in branching ratio calculation 
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FIG. 8: The Yoaiitz distribution for after all selection cuts. 
The dots show the data, the histogram the sum of the sig- 
nal and background Monte Carlo simulation, and the shaded 
histogram the scaled background Monte Carlo simulation. 



the kaon center-of-mass. The distributions of the Yoaiitz 
variable is shown in Figure[5J while the Zuaiitz variable is 
closely related to the m 3 4 distribution shown in Figure [5] 
The data used to determine the value of ay is listed in 



VII. SYSTEMATIC UNCERTAINTIES 

In general, the systematic uncertainties related to the 
branching ratio measurement are associated with either 
the acceptance calculation or the background estimate. 
The largest systematic error is due to the change in the 
acceptance as a function of the value of ay. This comes 
about mainly because the 11134 distribution depends upon 
the value of ay, and the acceptance varies across the 7)134 
region. We find that the ay dependence has the following 
form: B = (1.33 + 0.13 x ay) x 1(T 6 . We conservatively 
evaluate this systematic error by allowing ay to vary by 



±0.16. The next largest systematic uncertainties are re- 
lated to the acceptance ratio between the normalization 
and signal acceptances. To determine a systematic error 
for the acceptance, we compared the K L —> 7r°77 data 
and Monte Carlo simulation. We then reweighted the 
specific Monte Carlo distribution to match the same dis- 
tribution in data. We used this weight factor to calculate 
a new acceptance and used the difference to assign a sys- 
tematic uncertainty. We also examined the K L — > 7r°7r° 
decays and found similar results when reweighting the 
Monte Carlo simulation to match the data. Since the 
K L — ► 7T°7r and K L — > 7r°77 samples have nearly identi- 
cal acceptances, this gives us confidence in our estimate 
in the systematic effects. The variables that had the 
largest effect on the acceptance were the kaon energy 
and the photon veto response. We assigned a systematic 
uncertainty of 1.16% due to the acceptance. 

The ability of our Monte Carlo simulations to repro- 
duce the Kl — ► 7r°7r 7r background also contributes to 
the systematic uncertainty. To estimate the effects from 
our knowledge of the background, we looked at all events 
before applying the shape \ 2 cu t- This sample is dom- 
inated by Kl — + 7r 7r°7r° events. We then reweighted 
the background Monte Carlo sample to match the data 
in a particular parameter. The change in acceptance 
multiplied by the background fraction was taken to be 
the systematic uncertainty from a specific variable. In 
particular, we assigned the following systematic uncer- 
tainties due to our simulation of the background: the 
photon shape \ 2 (1-07%), the drift chamber simulation 
(0.92%), the photon veto simulation (0.90%), the kaon 
energy shape (0.69%), and the Kl decay vertex distribu- 
tion (0.38%). 

In addition to the acceptance calculation and the back- 
ground determination, a few other effects contribute to 
our systematic uncertainty including the Monte Carlo 
statistics, the background normalization and the mea- 
sured Kl — > 7r°7r° and 7r° — > 77 branching fractions. For 
this analysis, we generated nearly 10 11 Kl — > tt ^ ^ de- 
cays, more than twice the background statistics. These 
statistics contribute 1.0% to the systematic uncertainty. 
To determine the normalization of the Kl — * 7r 7r°7r° 
Monte Carlo sample, we first scaled the Kl — * 7r°7r° 
Monte Carlo sample to the observed 7r° peak in the 
TO34 mass distribution in the data. We then normal- 
ized the Kl — > 7r 7r°7r° sample relative to the number 
of Kl ~^ 7r°7r° events by the ratio of the branching ra- 
tios and the number of generated events. To assign a 
background normalization systematic error, we scaled the 
Kl — * 7r 7r°7r° background events directly to the shape 
X 2 distribution and compared the difference between the 
two methods. This contribututed 0.90% to the total sys- 
tematic uncertainty. Finally, we assigned a 0.5% sys- 
tematic uncertaintly due to the error on the measured 
Kl — > 7r°7r° branching ratio. All of the systematic effects 
are listed in Table HT1 with a total systematic uncertainty 
on the Kl —> 7r°77 branching ratio of 3.0%. 

To determine the systematic uncertainty in our ay 
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Type Source 


Uncertainty (%) 


Acceptance ay dependence 


1.50 


MC acceptance ratio 


1.16 


Background Photon shape \ 2 


1.07 


Drift chamber hits 


0.92 


Photon vetoes 


0.90 


Kaon energy 


0.69 


Decay Vertex 


0.38 


General MC statistics 


1.00 


Background normalization 


0.90 


Kl — * 2tv° branching ratio 


0.50 


Total 


3.0 



TABLE II: Branching ratio systematic uncertainties 



measurement, we varied the position of the selection cuts 
and looked for any non-statistical change in the value of 
ay. We also varied the level of Kl — ► 7r 7r 7r background 
according to the methods described above. The major 
systematic uncertainties associated with the determina- 
tion of ay are listed in Tabic IIIII The main sources of 
systematic error result from the uncertainty of the back- 
ground estimations. The total systematic uncertainty as- 
sociated with the ay measurement is 0.07. 



Source 


Uncertainty 


Z vertex cut 


0.05 


Photon veto cut 


0.04 


3-7T normalization 


0.03 


Photon shape x 2 


0.01 


Total 


0.07 



TABLE III: ay fitting systematic uncertainty. 



VIII. FINAL RESULTS AND CONCLUSIONS 

To obtain the final branching ratio result, we used 
the weighted average of the 1996-1997 and 1999 num- 
bers based upon the statistical errors of the two results. 
The systematic studies were done on the combined 1997 
and 1999 analyses to take into account any correlations. 
Including the uncertainties due to the systematic effects, 
we find the Kl — * 7r°77 branching ratio to be 

B{K L -> 7^77) = 

(1.29 ± 0.03 sfat ± 0.05 syst ) x 10~ 6 . (3) 

This result is a significant improvement over the previous 
KTeV result, and supersedes that result. The differences 
between the current and previous results is discussed in 
Section [Kl 

Our value of ay was obtained using the fitting method 
described above. The x 2 for the fit is 56.6 for 59 degrees 
of freedom. Including the systematic error, we find 

ay = -0.31±0.05 stat ±0.07 syst . (4) 



The total error from our determination of ay is slightly 
larger than the NA48 result, however, it is compatible 
with their value. 

The branching ratio result is consistent with the lat- 
est 0(p 6 ) ChPT results. Our value of ay suggests that 
the CP conserving amplitude in Kl — > n°l + l~ should 
be less than 1 x 10 -12 compared to the expected total 
branching ratio of ~ 3 x 10 _11 [10]. Therefore this de- 
cay should be dominated by CP violating terms. Future 
searches for Kl ^ n°e + e~ and Kl — > 7r°/i + ^~ would be 
of great interest since many models of new physics would 
signficantly alter these branching ratios. 



IX. APPENDIX A 

Compared to the previous KTeV Kl — * 7r°77 branch- 
ing ratio value, our new result is significantly lower. The 
main difference between the two analyses arises from our 
simulation of the transverse photon shower shape. Our 
previous analysis used the GEANT-based shower library 
while our current analysis utilizes the data-based shower 
library. As shown in Figure [31 the data-based shower 
library shows significant improvement over the GEANT- 
bascd shower library. Utilizing this new shower library 
changes our estimate of the Kl — > 7r 7r 7r° background, 
increasing the background by a factor of approximately 
two. The increase in background occurs because the 
Kl — ► 7r 7r°7r° background peaks in the region of small 
shape \ 2 as shown in Figure In our previous result 
the Kl — > 7r 7r 7r° background shape dropped in the sig- 
nal region. Accordingly, the background estimate uti- 
lizing the GEANT-based shower library underestimated 
the Kl — > 7r 7r°7r° background. 

In the previous result the systematic error for the mis- 
match in the photon shape x 2 scaled with the size of the 
estimated background. However, since the background 
shape was incorrectly modeled, our estimate of the sys- 
tematic error also was underestimated. Studies of other 
variables sensitive to the Kl — > 7r°7r 7r° background ar- 
rive at a similar estimate for the background level in the 
current analysis. 



Acknowledgments 

We gratefully acknowledge the support and effort of 
the Fermilab staff and the technical staffs of the par- 
ticipating institutions for their vital contributions. This 
work was supported in part by the U.S. Department of 
Energy, The National Science Foundation, The Ministry 
of Education and Science of Japan, Fundao de Amparo 
a Pesquisa do Estado de Sao Paulo-FAPESP, Conselho 
Nacional de Desenvolvimento Cientifico e Tecnologico- 
CNPq and CAPES-Ministerio Educao. 



8 



[1] G.D. Barr et al, Phys. Lett. B242, 523 (1990). 
[2] V. Papadimitriou et al, Phys. Rev. D44 R573 (1991). 
[3] G.D. Barr et al, Phys. Lett. B284, 440 (1992). 
[4] G. Ecker, A. Pich and E. de Rafael, Phys. Lett. B189, 
363 (1987). 

[5] G. D'Ambrosio and J. Portoles, Nucl. Phys. B492, 417 
(1997). 

[6] F. Gabbiani and G. Valencia, Phys. Rev D66, 074006 
(2002). 

[7] T.N. Truong, Phys. Lett B313 221 (1993). 

[8] A. Alavi-Harati et al. Phys. Rev. Lett. 93 021805 (2004). 

[9] A. Alavi-Harati et al Phys. Rev. Lett. 84 5279 (2000). 
[10] G. Buchalla, G. DAmbrosio and G. Isidori, Nucl. Phys 

B672 387 (2003). 
[11] F. Mescia, C. Smith and S. Trine, JHEP 0608, 88 (2006). 



[12] J.R. Batley et al, Phys. Lett B576 43 (2003). 

[13] A. Alavi-Harati et al, Phys. Rev. Lett. 83, 917 (1999). 

[14] A. Lai et al, Phys. lett B536, 229 (2002). 

[15] A. Alavi-Harati et al., Phys. Rev D67, 012005 (2003). 

[16] A. Alavi-Harati et al, Phys. Rev D67, 012005 (2003). 

[17] C. Bown et al, Nucl. Instr. Meth. A369, 248, (1996). 

[18] R. Brun and F. Carminati, CERN Program Library Long 
Writeup W5013 (unpublished). 

[19] T. Alexopoulos et al, Phys. Rev D70, 092006 (2004). 

[20] W.-M. Yao et al, J. Phys G33, 1 (2006). 

[21] Data used to determine both the Kt — > 7r°77 
branching ratio and ay will be uploaded as an 
EPAPS document. For more information on EPAPS, see 
http://www.aip.org/pubservs/epaps.html 



